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General Process for Developing and Implementing 
a Regulated Biochemical Reaction Network Model 
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Mathematically describe the regulatory controls 
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Formulate problem by specifying parameters and initial 
reaction constraints to the model 
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Example Biochemical Reaction Network 



B) Example Regulatory Structure and Requirements 



© S .0 



C) Simulated Reaction Network (without regulation) 



D) Simulated Reaction Network (with regulation) 
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Metabolic Model 



subj prod2 



Regulatory Model 



substrate, 
substratCj 
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•Activity constraints set for rxa^^^^^ 

(lower bound = 0, upper bound = INF or #) 

•Inactivity constraints for rxn^toj^j, 
(lower bound = 0, upper bound = 0) 



' prodj 



product, 
productj 



Integration of Stoichiometric model and 
Logical model achieved through regulatory 
restraints (logic values of reaction 
processes) which are used to refine 
appropriate reaction constraints in the 
model: 

If naiLogic^ 1 ths" "se Activity constraints 
'^Logic^ 0 then use Inactivity constraints 



Logic functions 

a, = (activator/inhibtor) • TF 

32=1 

Cj = TF* • pr, • genel • gene2 
^2 = pi's • geneS 

1, = Mg<.„^, 

Ij = Mg,„,2 

Pl ~ Pgaiel * Pgaie2 * PgeneS 

■^Logic^ Protein • Cofactor • Substrate, • Substrate2 

Time delays can be specified for the switching 
of each memorization variable after a 
triggering change in the associated function 
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A Time-Dependent Implementation of a Regulated 
Biochemical Reaction Network Model 
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Formulate problem by specifying the intiat time conditions 
and parameters 
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Refine the constraints on the biochemical reactions based 
on the status of their regulation 



p 



Use the flux distribution to detemnine the new set of intial 
conditions for the next time point 
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Fonnulate problem by specifying the conditions a 
parameters for the next time point 
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Solve for a steady-state flux distribution in the reaction 
network subject to the constraints 



END ^ 



Figure 6 



300 

Process for Developing Genome Scale Regulated 
Models of Metabolism 
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Obtain the DNA sequence of the organism strain 



Detennnine the location of all the open reading frames in the 
genome 



Determine genes involved in cellular metabolism and/or 
metabolic regubtion 
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Gather biochemical infomfiation on all reactions performed 
by the genes^ene products of the metabolic genoVpe 



Gather additional biochemical information on any reactions 
known to occur biochemically but not present in metabolic 
genotype 



Group metabolic and regulatory genes into transcr'ptional 
un'rts based on sequence anal^is 



Locate regulatory binding sites on transcriptional units by 
sequence ortnctional hortalogy 
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Specific default reaction constraint conditions based on 
themiodynamics and reaction kinetics where available 



Define the organism specific regulatory structure using 
Boolean logic representation 



Detemiine organism specific biochemical and physiological 

data on biomass composition, uptake rates and 
maintenance requirements 



Formulate the general linear programming problem 
combined VMth regulatory constraints representing the 
organism specific model 
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REACTION 




NAME 


REGULATION 


Metabolic Reactions 








-1 A-1 ATP+1 B 




R1 




-1 B +2 ATP +2 NADH + 


1 C 


R2a 


IF NOT(RPb) 


-1 C -2 ATP -2 NADH + 


B 


R2b 




-1 B +1 F 




R3 




-1 C+1 G 




R4 




-1 G + 0.8 C +2 NADH 




R5a 


IF NOT (RPo2) 


-1 G + 0.8 C +2 NADH 




R5b 


IF RPo2 


-1 C +2 ATP +3 D 




R6 




-1 C -4 NADH +3 E 




R7 


IF NOT (RPb) 


-16-1 ATP -2 NADH + 


H 


R8a 


IF NOT (RPh) 


+1 G +1 ATP +2 NADH 


1 H 


R8b 




-1 NADH -1 02 +1 ATP 




Rres 


IF NOT (RPo2) 


Transport Processes 








-1 Cailxinl +1 A 




Tc1 




-1 Cartxin2 +1 A 




Tc2 


IF NOT(RPcl) 


-1 Fext +1 F 




Tf 




-1 D +1 Dext 




Td 




-1 E +1 Eext 




Te 




-1 Next +1 H 




Th 




-1 Oxygen +1 02 




To2 




Maintenance/Growth Pr 


xesses 






-1 C -1 F -1 H -10 ATP +1 Biomass 


Growth 




Regulatory Proteins 












RPo2 


IF NOT(Oxygen) 
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IF Carboni 
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IFTh 
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